package scu.maqiang.tiangong;

import scu.maqiang.fes.BVPType;
import scu.maqiang.mesh.Mesh3T4;
import scu.maqiang.mesh.MeshFormat;
import scu.maqiang.mesh.ScalarFunc;
import scu.maqiang.possion.Possion3DT41;

public class ZKEngineMeshTest {

	public static void main(String[] args) {
		String dir = "TianGong\\ZKEngine\\";
		String meshFile = dir + "ZKEngine.mesh";
		Mesh3T4 mesh = new Mesh3T4().readMesh(meshFile, MeshFormat.mesh);
		
        ScalarFunc cf = (xyz, label, param) -> 1.0;
        ScalarFunc sourceFunc = (xy, llabel, param) -> 10.0;
        Possion3DT41 p2dt31 = new Possion3DT41(mesh);
        p2dt31.DirichletProblem(BVPType.COMMON, cf, null, sourceFunc, null,
        		                null, null, null, dir + "ZKEnginePossionSolution.dat");
	}

}
